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Abstract 

We calculate non-perturbative renormalization factors at hadronic scale for AS = 2 four-quark opera- 
tors in quenched domain-wall QCD using the Schrodinger functional method. Combining them with the 
non-perturbative renormalization group running by the Alpha collaboration, our result yields the fully 
non-perturbative renormalization factor, which converts the lattice bare Bk to the renormalization group 
invariant (RGI) Bk- Applying this to the bare Bk previously obtained by the CP-PACS collaboration at 



a. 



t 2,3,4 GeV, we obtain B K = 0.782(5)(7) (equivalent to B% S (NDR, 2GeV) = 0.565(4)(5) by 2-loop 
running) in the continuum limit, where the first error is statistical and the second is systematic due to the 
continuum extrapolation. Except the quenching error, the total error we have achieved is less than 2%, which 
is much smaller than the previous ones. Taking the same procedure, we obtain mj^ 1 = 5.613(66) MeV and 
to rgi = 147.1(17) MeV (equivalent to m^f(2GeV) = 4.026(48) MeV and mp(2GeV) = 105.6(12) MeV by 
4- loop running) in the continuum limit. 
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I. INTRODUCTION 



In the standard model, the dimension-six four-quark operator, 

LL = s^(l- l5 )d-s llM (l-^)d, (LI) 

of the low-energy effective Hamiltonian induces the K° — K° mixing, and the estimation of its 
hadronic matrix element (K°\Oll\K } is required to extract Cabibbo-Kobayashi-Maskawa matrix 
elements from the experimental value of the indirect CP violation parameter ex- The hadronic 
matrix element is parametrized by the kaon B parameter Bk, defined by 

B _ (A- |^(l-75)d-^(l-75)^°) n2) 
K (8/3)(i?0|s7 M 75d|0)(0|s 7/t 75d|^ > ' K ' 1 

and lattice QCD can provide the first principle calculation of it. In the past decades much effort 
have been devoted to the estimation of Bk by employing various quark and gauge actions [H, 0, 0, 
0, @, 0, S, S E3, 11, 12, Recently it is recognized that an essential step toward the precise 



determination of Bk is to control the systematic error associated with the renormalization, and for 
the precision now required, the non-perturbative renormalization seems necessary 0, EH. 13], 



Among several non-perturbative schemes on the lattice the Schrodinger functional (SF) scheme 
[13 . Ea . 0, E3] has an advantage that systematic errors can be unambiguously controlled: A 
unique renormalization scale is introduced through the box size to reduce the lattice artifact and a 
large range of the renormalization scale can be covered by the step scaling function (SSF) technique. 

A few years ago the CP-PACS collaboration has calculated Bk using the quenched domain-wall 
QCD (DWQCD) with the Iwasaki gauge action [g], and a good scaling behavior with small statis- 
tical errors has been observed. Systematic errors associated with the perturbative renormalization 
factor at one loop, however, can not be precisely estimated. A main purpose of this paper is to 
remove this uncertainty of the renormalization factor, by evaluating it non-perturbatively. 

We adopt the SF scheme to control systematic uncertainties due to the finite lattice spacing. 
In the SF schemes, the renormalization factor Zs K (go), which convert the bare Bk to the renor- 
malization group invariant (RGI) Bk, is decomposed into three steps as 



(1.3) 

at a given bare coupling. 

The first one is the renornialization factor at the hadronic scale /-/-min; 

which is given by 



Z 2 A (9o) 



Bjc (SO , a/^min) = ^PTT~\ ' ( L4 ) 



where Zyv+AA and Za are the renormalization factors for the parity even part of Oll an d for 
the axial vector current, respectively. Here a/x m i n <C 1 should always be satisfied to keep the 
lattice artifact small enough for the reliable continuum extrapolation. This factor depends on both 
renormalization scheme and lattice regularization. Multiplying it by the lattice bare operator, the 
regularization dependence is removed and only the scheme dependence remains. 

-^FA+AV (/•'max, Mmix) represents non-perturbative RG running for the parity odd part of Oll, 
from the low energy scale //mm to the high energy scale /i max = 2 7 /U m i n where perturbation theory 
can be safely applied. Among three steps this part requires the most extensive calculation. Since 
this factor does not depend on a specific lattice regularization after the continuum extrapolation, 
we can employ Zy\ +Av (ij, 

max, /^min) evaluated previously by the Alpha collaboration with the 
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improved Wilson fermion action [18|], instead of calculating it by ourselves. Note that the renormal- 
ization factors for the parity even and the parity odd parts agree after the continuum extrapolation, 
thanks to the chiral symmetry. 

The last factor Zy^ +Av (oo, u max ) is the RG evolution from the high energy scale /U max to 
infinity, which absorbs the scale dependence to give the RGI operator. Since we are already deep 
in the perturbative reg ion at fimax, we can evaluate this factor perturbatively, using the two loop 



calculation in Ref. 121]. Note that the scheme dependence is also removed at this stage and the 
RGI operator becomes scheme independent. 

Our target in this paper is the calculation of the first factor Z^^(go, a/i m in)- I n order to further 
reduce the computational cost^we use a relation that Zy = Z A implied by the chiral symmetry of 



the DWQCD in SF scheme [20], together with another one that Zyv+AA = Zva+av^ which will 
be checked numerically in this paper. Therefore, throughout this paper, we adopt the following 
definition, 

This paper is organized as follows. In section[II]we introduce the SF renormalization scheme and 
RGI operator for Bk following the Alpha collaboration. Numerical simulation details are described 
in section IIII1 In section IIVI we present our main results for the non-perturbatively renormalized 
RGI Bk, and we discuss its continuum extrapolation. We have also made several numerical checks 
of our formulation. Section [V] is devoted to the non-perturbative renormalization of light quark 
masses. Our conclusion and discussion are given in section IVTl 

II. SCHRODINGER FUNCTIONAL SCHEME AND RGI OPERATOR 
A. Renormalization group invariant operator 

A bare n-point correlation function on the lattice, 

G (xi,--- ,x n ;a;g ,m ) = (Oi(xi) • • • O n {x n )) , (II.6) 
is multiplicatively renormalized in the mass independent scheme as 

Gr{x x ,--- ,x n ;n;g R (n),m R (fi)) = ^JJ Z 0i (go, o/i)^ G (xi,--- , x n ; a; g Q , m ), (II.7) 

where gu and ttir are the gauge coupling and the quark mass respectively, while corresponding 
bare quantities have the subscript 0. 

The RG equation for the n-point function reads 

(d d d n \ 



where 



(II. 



P{9R) = ^~d^T' t{9r) = mM ^ ' (IL9) 

10A9R) = lim ~^— ( r/x t. • (11.10) 

a^o Z 0l {go,afi) a ii 
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From the RG equation, the finite scale evolution of G R from \x to //' is calculated as 

G R (xi,--- ,x n ;n';g R (p'),m R (fi')) = JJ [/*(//, /z) G R {xi,--- ,x n ; fi;g R (fji),m R (pi)), (11.11) 



where 

Ui{p,p) = exp / dg , r = lim — 11.12) 

is the scale evolution for each operator in the continuum limit. Using this factor we can define the 
RGI operator 0(x) as 

(0) 



-»(-/„ *(w-t?jj « ( ^ <II13) 

where bo and 7q are given by 

/?(<?) = -b 09 3 - hg 5 - b 2 g 7 + ■■■ (II. 14) 

7o(g) = -7SV - 7 gV - 7 gV + • • • , (H.15) 

and O r (x; fj,) is the renormalized operator at some scale /i. 

As mentioned in the introduction, the evaluation of the RGI operator in the SF scheme is decom- 
posed into three steps. The lattice bare operator is renormalized at scale /u m i n non-perturbatively 
with the first factor ZQ P (go,ajj, m - m ). The scale evolution from /x m i n to /i max is given by the second 
one, 

■^0 P (/ i max5 A'min) = U(/J, nmx , /i m in) > (11.16) 

which can be evaluated non-perturbatively using the step scaling function. The last factor is the 
running from /i max to infinity, which can be calculated safely by the perturbative expansion as 



„t°) 

/„2 (,. \\ — 



B. Schrodinger functional scheme 

In the SF method, the renormalization scheme is specified by the choice of the correlation 



function in the finite box. Since we rely on the result by the Alpha collaboration 18[] for the RG 
running from /i m in to fi max , the same correlation function must be taken as the renormalization 
scheme for our definition of Z^(go, fi m i n ). 

We here consider The following form of the correlation function, 

F± aTbTc (x ) = ^(0 21 [T A }O i5 [T B }0$ A+AV (x)0' 53 [r c ]), (11.18) 
where subscripts 1 ~ 5 represent quark flavours, and 

°VA+AV = 2 [i^^M (^37^75^4) + (^17^75*02) (^37^4) 

± ((^17^4) (^37/^75^2) + (^17/^75^4) (^37^2))) (11.19) 
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is the parity odd four quark operator made of four different flavors. Boundary operators Oij and 
0' i: j are given in terms of boundary fields C and (f [17| as 

O^T] = a 6 ^ C,(.r)\X j{ <f). OijiT] = a 6 ^ ( IL2 0) 

xy xy 

Due to the SF boundary condition for fermion fields the boundary operator should be parity odd 
and we then have two independent choices, r = 75 and r = 7^ (k = 1,2,3). For the correlation 
function to be totally parity-even we need at least three boundary operators as in ()II.18j) . 
The Alpha collaboration has adopted five independent choices for the correlation function 

Ff(x ) = ^ 5i75i75 (x ), 



j,k,i=i 

F ?(*o) = ^E^, M (^o), (11-21) 
fc=i 

1 3 

fe=i 
1 3 



3 ,=1 

To remove logarithmic divergences of boundary fields £'s from these correlation functions, one can 
consider the following 9 ratio of the correlation functions, 

hf(xo) = ^S^l, i = l,...,5, 
fi 

h£i*o) = (H.22) 

hf +i (x ) = 4 = 3,4,5, 

A fa 



where 

1 1 3 

h = -^{o' 12 h 5 ]o 21 [ lB }), h = --^Y,( 'M°M), (H.23) 

fc=i 

are the boundary-boundary correlation functions. Each ratio, distinguished by the label s = 
1, • • • ,9, gives a different renormalization scheme. Among these 9 choices, /i^ =137 (xo) and 
^7=2 4 5 6 8 s( x o) define good schemes [l^], whose scaling violations are perturbatively shown to 
be small. 

The renormalization factor we need in this study is defined by 

ZvA + AV;s(9a,»)hf(xo = L/2;g ) = hf^(x = L/2), (11.24) 

where s labels the scheme, and ht^ tree ^ is the correlation function at the tree level in the continuum 
theory. 
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According to Ref. [21[ , the renormalization factor for the local vector current is defined through 
the Ward-Takahashi identity as 

Zv = f v (xo = L/2y fv{Xo) = ? < O 'i2[75]Vb(x,x )O3i[75]> , (H.25) 

where V^(x) = ip 2 ( x )l^ V>3 (?) ■ 

III. NUMERICAL SIMULATION DETAILS 
A. Gauge action 



The theory is defined on an L 3 x T lattice of L = T = N L a = N T a [18|], with the periodic 
boundary condition in the spatial directions and the Dirichlet boundary condition in the temporal 
direction. The dynamical gauge variables are spatial links Uk(x) at xq = 1, . . . , Nt — 1 and temporal 
ones Uq(x) at Xq = 0, . . . , Nj- — 1. The Dirichlet boundary condition is imposed on the spatial link 
at xo = and Nt as 

U k (x,x = 0) = exp[aC fc ], U k (x,x = N T ) = exp[aC' k ], (111.26) 



where C k and C' k are anti-Hermitian diagonal matrices [1J, [15|] , which we set to zero in our simu- 
lation. 

We employ the renormalization group improved gauge action , 

S g iuo„[t/] = - I W P (gl) Yl Retr ( J - ^pi) + w R(9o) E Retr ( J " U ^ I ( IIL27 ) 



(Jo 



plaquette rectangle 



where U p \ represents the standard plaquette and U r tg an 1 x 2 six-link rectangle. The 0(a) lattice 
artifact due to temporal boundary is removed by setting weight factors Wp and Wr as 

2 \ J cocflfi'o) Set of temporal plaquettes that just touch one of the boundaries 
cq Otherwise 



w R ( 9 l) 



cic^((/q) Set of temporal rectangles that have exactly two links on a boundary 
ci Otherwise 



The coefficients cq and c\ are normalized such that cq + 8ci = 1. In this paper we take c\ = —0.331 
(the Iwasaki gauge action) . The boundary coefficients are expanded perturbatively as 

cocf (So 2 ) = co (l + cf (1) 5o 2 + 0(gt)) , (111.28) 



Cl cf (gl) = Cl ( 1 + cf {l) gl + 0(gt) ) . ( 111.29) 



Since only a single improvement condition 



c cf (1) + 4cicf (1) = 0.1518 (111.30) 



is available [231] . there exists no unique choice. Therefore, in this study, we adopt the condition 



A m that cf (1) = 2cf (1) 
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B. Fermion action 



In this paper we adopt the orbifolding construction of the SF formalism for the domain-wall 
fermion 25|, |2g, |27]] . Instead of folding the temporal direction as was discussed in Ref. 27j], we keep 
both positive and negative regions in the temporal direction, in order to implement the even-odd 
preconditioning in five dimensions. The gauge link in the negative region is defined to satisfy the 
time reflection symmetry as 

U k {x,x ) = U k (x,-x ), Uq(x,x ) = U^(x, -xq - 1). (111.31) 

We implement the Shamir's domain-wall fermion action 2^, 2^] on 2Nt X x N$ lattice, 

N T N 5 

Sdwf = ^2 X] ^2 ^( x i s ) D dw{{x,s;y,t)Tp(y,t), (111.32) 

x,y x ,y =-N T +l s,t=l 

where the temporal coordinates xq and yo run from — Nj- + 1 to Nt, while the fifth dimensional 
coordinates s and t from 1 to N§. For the orbifolding we set the anti-periodic boundary condition 
in the temporal direction, 

1p(x, Xq + 2Nt, s) = — 1p(x, Xq, s), 1p(x, Xq + 2iVy , s) = — 1p(x, Xq, s) . (III. 33) 



On the other hand, the periodic boundary condition with the phase = 1/2 in spatial directions [la. 



1p(x k + N L , Xq, s) = e l9 1p(x k , Xq, s), 1p(x k + N L , Xq, s) = e~ l9 1p(x k , Xq, s), (111.34) 

is imposed by replacing the spatial gauge link as Uk(x) —> e %e f NL Uk{x). We set the physical quark 
mass to be zero for the mass independent (massless) scheme. 
The physical quark field is defined in the standard manner as 

q{x) = (P L S s ,i + Pr5 SjNb )iP(x, s), q(x) = ^{x, s)(S s , Ns P L + S Sjl P R ), (111.35) 

with Pl/r = (1 ± 75)/2, and its propagator on the 2^-^ x iV"£ x ./V5 lattice is given by 

G quark (x, y) = (P L S S>1 + P R 5 S , N5 ) (-=}—) (S t>NB P L + 5 tA P R ) . (111.36) 

V^dwf / x ,y;s,t 

Imposing the orbifolding projection we get the physical quark propagator in the SF formalism as 

GgLrkfoi/) = 2(n_G quark n + ) ;C)2/J n ± = (111.37) 

where R is the time reflection operator: R XQ)yo q(x,yQ) = q(x,—xo). Due to the projection, the 
physical quark fields satisfy the proper homogeneous SF Dirichlet boundary condition at xq = 0, Nt 
such that 

P + q(x)\ Xo=0 = 0, P-q(x)\ Xo=NT =0, (HI.38) 
q(x)P.\ Xo= o = 0, q(x)P + \ XQ=NT = 0. (111.39) 

As usual, the boundary-bulk and boundary-boundary propagator are constructed in terms of the 
SF quark propagator pil.37[) [la ]. 
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C. Parameters 



The CP-PACS collaboration has calculated the lattice bare Bk in quenched DWQCD with the 
Iwasaki gauge action at the domain wall height M = 1.8 and the fifth dimensional length N§ = 16 
0]. In order to renormalize this Bk we have to take the same lattice formulation. The bare value 
of Bk, calculated at three lattice spacings (3 = 2.6, 2.9 and 3.2* (a -1 ~ 2, 3 and 4 GeV) in the 
previous simulation [6], is listed in tabled 

The renormalization scale at the low energy (hadronic scale) is introduced as l/yu m in = 2Lpv, 
where L max is defined through the renormalized coupling g 2 (l/L max ) = 3.480 in the SF scheme [181 ]. 
and I/maxAo = 0.749(18) [24| in the continuum limit. (/« m m = l/2£ m ax ~ 263 MeV for tq = 0.5 
fm). At (3 = 2.6, 2.9 and 3^2, N L which satisfies aN L = 2L max — 1.498ro can be estimated, using 
the interpolation formula [241 ] . 



In 1-1 = -2.193 - 1.344(/3 - 3) + 0.191(/3 - 3)'*, (111.40) 

valid at 2.456 < (3 < 3.53. To cover the resulting lattice sizes, N L = 7.60625, 11.7144, 17.4317, we 
take 7 lattice sizes, Nl = 6, • • • , 18, and using the formula (jIII.40jl again, we tune (3 so that the 
physical box size satisfies clNl — 2i max — 1.498r at each N L . 

Quenched gauge configurations are generated by the HMC algorithm. First 2000 trajectories 
are discarded for thermalization, and the correlation functions are calculated every 200 trajectories. 
By the jackknife analysis we found that each configuration separated by 200 trajectories is almost 
independent. A value of (3 and a number of configurations at each lattice size are listed in table ITU 



IV. NON-PERTURBATIVE RENORMALIZATION OF B K 
A. Extraction of renormalization factors 

The behavior of Zy given in (1TL251) is plotted function of time xo in Fig. [1] at Nl = 6 ~ 20. 
As a decreases (./Vi increases), Zy becomes flatter in xq. Typical behaviors of Zy A+Ay . s (|II24|) is 
given for the schemes s = 1 in Fig. [2] and the schemes s = 2 in Fig. [3j and Zy A+Av . s for s = 1 
in Fig. [Hand s = 2 in Fig. \5\ Both renormalization factors are almost xq independent for s = 1, 
while they strongly depend on xq for s = 2. 

Taking the value at xo = L/2[la|, we get renormalization factors, whose numerical value are 
listed in table IIIII Combining Zy A+Ay . s and Zy, we get the renormalization factors for Bk in (|I.5p . 
which is also listed in the table IIIII All errors in the table are evaluated by a single elimination 
jackknife procedure. 



B. Scaling behavior of the step scaling function at L max 

In this subsection, we discuss universality of the scale evolution function Zy A , ^(/Jmaxi Mmin)- 
More explicitly we calculate the SSF at the largest coupling u = g 2 (l/L max ) = 3.480 for four values 
of f3's (lattice spacings) to make the continuum extrapolation, and compare the result with that 
by the Alpha collaboration [3]. The SSF of the four-fermion operator is defined as a ratio of the 



* The data at = 3.2 is new and not published in Its numerical analysis is briefly given in appendix iBl 
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renormalization factors at two different box sizes: 



^VA+AV;s( U > a / L ) 



z 



VA+AV; 



,(S0,a/(2L)) 



7 ± 

VA^ 



AV;s(9o,a/L) 



m=0,g 2 (l/L)=u 



(IV.41) 



Since we have already calculated Z ± (go, aj (2L max )) in the previous subsection, we need to calculate 
Z ± (go,a/L max ) except Nl = 20 = 2L max /a. Number of configuration is fixed to 100 for all Nl- 
Values of Z ± (go, a/L max ) at Nl = 4,6,8, 10 and Z ± (go, a/(2L max )) at Nl = 20 are given in table 

ED 

In Fig.El we compare our SSF (filled circles) for , av-\( u i a /^ J ) (Left) and ^va+AV-s( u ^ a /^) 
(Right) with those by the Alpha collaboration calculated with the clover fermion (triangle up) and 
the Wilson fermion (triangle down) as a function of a/L, together with their combined continuum 
limit by star 18|]. We have surprisingly found that the scaling violation of our SSF is large and 
they seem to approach their continuum limits with oscillation. To check whether this oscillating 
behavior is caused by the 0(a) bulk chiral symmetry breaking effect of the DWQCD at finite N$ or 
not, we investigate the A% dependence of the SSF at u = 3.480. As is shown in Fig. comparisons 
between A% = 8 and A5 = 16 for HyA+AV;i( u i L/a = 4) and between A% = 32 and A% = 16 for 
T,VA+AV;i( u , L/a = 6) indicate no A% dependence within statistical errors. The O(o) bulk chiral 
symmetry breaking effect has nothing to do with the oscillating behavior. 

We then suspect that the bad scaling behavior is caused by the 0(a) boundary effect in the SF 
scheme of DWQCD. To confirm this, we calculate the tree level SSF on the lattice, 



J VA+AV 



ML) 



Z VA+AV;l( a / 2L ) 



lattice 
tree 



z 



VA+AV 



M L ) 



lattice 
tree 



(IV.42) 



where N$ — ► oo limit is already taken. At tree level, we have £ 



VA+AV; 



ML) 



J VA+AV;8 



(a/L), 

which of course approach to 1 in the continuum limit. In this calculation, we take the tadpole 
improved value M ta( j = 1.5 at /3 ~ 2.9 for the value of M instead of the tree level value M = 1.8, 
in order to take into account an additive shift of M caused by the quantum correction. We plot 
the scaling behavior of ^>vA+Av-i^ a /^ open circles in Fig. [HI which shows an oscillation similar 
to one in Fig. [6l On the other hand, if we take M is close to but smaller than unity, the scaling 
behavior is much improved without oscillation, as is shown by open triangles at M = 0.9 in Fig[8j 

The tree level analysis indicates that the scaling behavior can be improved by changing the 
domain-wall height M so that the tadpole improved value becomes close to unity. Motivated by 
this, we have recalculated the non-perturbative SSF at M = 1.4, which corresponds to M ta d — 1.0 
at the range of our (3. Results are given in table [V] and are plotted by open diamonds in Fig[9l It 
is clearly seen that the scaling behavior at M = 1.4 is much improved, so that the linear continuum 
extrapolation can be made using last three points. The value in the continuum limit (filled symbol) 
is consistent with the previous one by the Alpha collaboration (star). 

We explore a different method to improve the scaling behavior of the SSF, without performing 
new simulations at different value of M. A main idea is to cancel the oscillating behavior of the 
SSF by that at tree level, changing the renormalization condition from 01.24(1 to 



7 ± 

^VA+AV-^s 



(go,n)hf(xo = L/2;g ) = h 



±(tree) 
s (lattice) 



(x = L/2), 



(IV.43) 



where the tree level correlation function, evaluated at M = M ta d for corresponding (3, is used in 
the right-hand side. We call this method the tree level improvement. Results are given in table IVT1 
and IVIH and are plotted by open squares in the Fig [9] We find that the magnitude of oscillation 
is reduced, so that a linear continuum extrapolation using last three data becomes possible. The 
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value in the continuum limit is consistent with both one by the Alpha collaboration and one at 
M = 1.4. 

In addition to the SSF of V A + AV, we have also considered the SSF of Bk defined by 

V I IT\ Z B K (90,a/2L) (T\T AA\ 

Y< BK (u,a/L) = — — — (IV.44) 

Zb k {go, a/L) 

at L = L max . Results are plotted as a function of a/L in Fig. [10] for three "good schemes". In 
each figure, results at M = 1.8 with and without the tree level improvement are represented by 
open squares and open circles, respectively, while the result at M = 1.4 by open diamonds. The 
scaling behaviors are reasonably well with the tree level improvement or at M = 1.4. Even without 
improvement, the oscillation is not so large. Linear extrapolations with three data at finest lattice 
spacings give consistent results among all three cases. The large oscillating behavior seems to be 
partly canceled between Zy A+AV and Zy in Zb k ■ 

We finally perform combined linear fits of the M = 1.4 data and the tree level improved data 
using the finest three lattice spacings. Values in the continuum limit of all SSF are given in table 

Em 



C. Renormalization of Bk 

In this subsection we evaluate the renormalization factor 2b k - s (9q) which convert the lattice 
bare -Br-(<?o) of DWQCD to the RGI Bk- As suggested in the previous subsection, we here employ 
the renormalization factor obtained with the tree level improved condition, hoping that this also 
improves the scaling behavior of Bk- Combining our renormalization factor Z|j P . s (<7o, A'min) in 
table ED with the RG running factor Zy A _^ A y(oc, /^max)-^y5+Ay(/^ max ' 

A*min) given by the Alpha 
collaboration [l~8l |. we obtain the renormalization factor Zb k ;s(9o) at each (3 in table iDCl 

In order to obtain the renormalization factors at j3 = 2.6, 2.9 and 3.2, we interpolate the result 
at each scheme s by the polynomial, 

Z Bk -M = a s + 6,09 - 3) + c,(P - 3) 2 , (IV.45) 

which is shown in the left panel of Fig. [TT] for schemes 1, 3, 7, together with interpolated values at 
three /3's by solid symbols. As is shown in the figure the polynomial interpolation works very well 
with small x 2 /dof ~ 0.1. Since the renormalization factor Zb k - s should not depend on schemes, 
the discrepancy between three schemes is considered to be the lattice artifact and therefore it 
should disappear at high /?, as seen in the figure. The renormalization factors with and without 
the tree level improvement are compared in the right panel of Fig. [TT1 Two renormalization factors 
are consistent within statistical errors, and agree completely at high (3 as expected. 

Multiplying the renormalization factor to the bare BK(go), we obtain the RGI Bk in table IXl 
The scaling behavior of Bk is shown in the Fig. [12] for s = 1,3,7, as a function of my a. Note 
that the scaling behavior of Bk with other (bad) schemes is indeed bad, therefore we do not use 
them in our analysis. Since the scaling violations are small, we have made the constant continuum 
extrapolation using the last two data points. We arrive at 

B K = 0.783(9) for Scheme 1, (IV.46) 
B K = 0.776(10) for Scheme 3, (IV.47) 
B K = 0.786(9) for Scheme 7. (IV.48) 

Since the values in the three schemes agree within errors in the continuum limit, we have made 
the combined constant fit for all three schemes, which gives 

B K = 0.782(5) (x 2 /doi = 0.87). (IV.49) 
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To estimate the ambiguity of the continuum extrapolation, we have also made the combined linear 
extrapolation using all 9 data points, and we obtain 



B K = 0.789(14) (x 2 /dof = 0.16). 
Now the final result we obtain leads 

B K = 0.782(5)(7) 



(IV.50) 



(IV.51) 



where the central value and the first error are taken from the combined constant fit, while the 
systematic error, given in the second, is estimated by the difference between the constant and 
the linear fits. Our result is consistent with previous results nonperturbatively renormalized by 



DWQCD [ll| with the RI/MOM scheme (B K = 0.786(31)) and by tmQCD [lj, \l$ with the SF 

scheme (B K = 0.735(71)). 

For the latter convenience, we convert Bk to the renormalized Bk in MS scheme with the naive 
dimensional regularization (NDR) at a scale fj, = 2 GeV. The renormalized operator in MS scheme 
is obtained by inverting the definition of the RGI operator as 



Air 



exp 




(o)' 
To 

bog 



O(x), 



(IV.52) 



where /%jg(ff) and 7q S (#) are renormalization group functions in MS scheme, which are estimated 
at four loops for /%js(<?) [20] and at two loops for 7q S (<7) 3l|. The gauge coupling g^(fi) in MS 
scheme is given in terms of Aj^g as 



A Ms = MOo0ns(^)) 2 "° exp 



1 






exp 



s: 



r(M) 



dg 



1 



+ 



1 



(IV.53) 



We adopt a value A^jg = 0.586(48)/Vo in Ref. [32], and take ro = 0.5 fm to set a scale. Multiplying 
Zb k -, s (9o) with the factor in (|IV.52p . we obtain the renormalization factor Zg^(go, 2GeV), which 
is listed in table IXI1 The scaling behavior of B^ s is given in the Fig. [T3j Note that our non- 



perturbative result differs from at (3 = 2.9 but agrees with at (5 = 3.2 with the previous result 
of the DWQCD with perturbative renormalization 33j, [34j, [35J] . The continuum extrapolation has 
been made as before, and we obtain the final result, 



_B^ S (NDR, 2GeV) = 0.565(4)(5). 



(IV.54) 



D. Chiral symmetry breaking effect 

As a further check, we investigate whether the assumption that Zyy+AA = %va+av holds or 
not in our DWQCD. We give only a result here and a detail of the SF formalism relevant in the 
analysis can be found in appendix A. 

If the chiral symmetry were exact, under chiral rotation of the first flavor 

qi ->• qx = *759i, Ci->Ci = -*75Ci> Ci -> Ci = -i7sCi, (IV.55) 
we could have the chiral Ward-Takahashi identity, 

(O va+av O[C})s = {Ovv+aaO{C]) s , (IV.56) 
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where 0[Q is the chirally rotated boundary operator of pi,20|) . A subscript S represents an action 
under which the expectation value is evaluated. Note that the boundary fields are also rotated, 
which satisfies opposite SF boundary condition to (1III.39I) . Therefore both actions are identical 
in the bulk but have opposite temporal boundary conditions. From this WT identity we obtain 
Zw+AA = Zva+AV- 

Unfortunately the domain-wall fermion action has a non-invariant part under the chiral rotation 

as 



S dwi - S dwl + Y, (IV. 57) 

where Y = ipXip is the bulk chiral symmetry violating term at the middle of the fifth dimension. 
Therefore the WT identity becomes 

(O va+av O[C})s = (O vv+aa 6[(})s + y + (O vv+AA 6[(})s. (IV.58) 

A possible chiral symmetry violation comes from the contribution of Y, which is expected to be 
suppressed exponentially in N$. We estimate the violating effect by comparing (Oy A+A yO[Q\) s 
with (0 V v+aaO[(]) s directly. 

We evaluate the renormalization factor Zyy +AA (go, /i m i n ) using the chirally rotated boundary 
condition and correlation functions of (|II.21|) and (111. 231) with the statistics of 100 configurations. 
The results are listed in table IXIII for the unimproved renormalization condition, and the time 
dependences of Zy A+Ay and Zy V+AA are shown in Fig. [HI and [15] for schemes 1 and 8, respectively. 
We observe good agreements between them at all L/a, and a similar results are obtained at other 
schemes. This investigation concludes that the relation Zyy +AA = Zy A+A y holds within statistical 
errors in our simulations. 



V. NON-PERTURBATIVE RENORMALIZATION OF QUARK MASSES 

From the relation derived from the axial- vector Ward-Takahashi identity in DWQCD^j that 

Z m = (V.59) 

we obtain the renormalization factor of quark masses from Z"p, which can easily be extracted as 
a by-product of the calculation in the previous sections. In this section, we report our results for 
the nonperturbative renormalization of quark masses. 



A. Renormalization group invariant quark mass 



The renormalization group invariant quark mass is defined by 

M = m(ji) (26 3 2 Gu)r^ exp (- dg (J^ - V (V.60) 

where m(fj,) is a renormalized mass in some scheme at scale [i. We evaluate the renormalization 
factor Zm, which converts the bare quark mass on the lattice in DWQCD to the RGI quark mass. 
A strategy to derive the renormalization factor is the same as that for Bk, and we write 

Zm{90) = Z^ioO, /X max )Z^ P (^ max , /X min )Z"^ P (5o, OMmin)- (V.61) 
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The first two factors have already been calculated by the Alpha collaboration as 3^, 37] 

< T (oo, aw)^ P (/W, /; min ) = _ M = 1.157(12) (V.62) 

at the same scale fj, m - m = l/(2L max ) as B^. As in the case for Bx what we need to calculate is the 
third factor, 

Z m (go, Q/Xmin) = ry , r- (V.63) 

The Alpha collaboration adopted the definition for the renormalization factor of the pseu- 
doscalar density such that 

Zp{l/L) = lA^Tj2Y f^) = -l( P S^o,S)O a ), fl = --L(O' a O a ), (V.64) 

where 

P a (x) = q(x)^r a q(x), O a = a 6 Y^ C»7 5 r a C(y) (V.65) 

xy 

are bulk and boundary pseudoscalar densities. 

Parameters for numerical simulations are same as those in section IIII CI A typical behavior of 
Zp as a function of time is shown in Fig. [TBI Values of Zp and Z m at xq = L/2 are listed in table 
IIIH with errors evaluated by a single elimination jackknife procedure. 

B. Scaling behavior of the step scaling function at L max 



We again study the scaling behavior of the SSF, 

Zp{go,a/2L) 



£p(it, a/L) 



Zp(go,a/L) 



(V.66) 

m=0,g 2 (l/L)=u 



at L = L max , which can be calculated from data in table ITO1 and HVl As seen in Fig ll7l the scaling 
violation in this case is also large and it seems to approach the continuum limit with oscillation. 

As before, we try to improve the scaling behavior by either taking M = 1.4 or using the tree 
level improvement with the renormalization condition that 



/p(xo = L/2) (tree) /p(xq = £/2)jg t t e) 

V ^^(liittice) 



(V.67) 



where (Zp)^ v ^) s = — — — -. The SSF obtained from the ratio of Zp in tables |V] [VlJ 

lattice; ^ _ M tad ) 

and IVIll is plotted in Fig. [18] for data at M = 1.4(open diamonds) and for the tree level improve- 
ment (open squares). In both cases the scaling behaviors are improved and the linear continuum 
extrapolation using the finest three data in each case are consistent with the value by the Alpha 
collaboration. A combined linear fit to both data gives CF m (u) = 0.853(13) with ^ 2 /dof = 0.71. 
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C. Renormalization of quark masses 



We employ the tree level improved condition (|V.67p for the renormalization of quark masses. 
Multiplying Z^ l p (gQ, a// m j n ) in table [VI] with the RG running factor <|V.62|) . we obtain the renor- 
malization factor Zm in table IXIIIl which is plotted in Fig. [19] by open triangles, together with 
data in table IXIVI at (3 = 2.6, 2.9 and 3.2 (filled symbols) by the quadratic interpolation. Data 
without the improvement are also shown in the figure. A discrepancy between the two are clearly 
observed at low (3. 

We convert the bare masses in table [I] to the RGI light quark masses, which is listed in table 
IXIVI Here m u d is the up and down averaged quark mass determined by it, while m s (K) or m s ((f>) 
is the strange quark mass by K or </>, respectively, and — m res is the residual mass of the DWQCD 
at which the pion mass vanishes. Following the previous paper Q], we adopt m q + m res for our 
definition of quark masses. Lattice spacing is given with p meson input. 

Since the scaling behavior of the RGI quark masses is reasonably good as shown in Fig. [20] we 
take the constant continuum extrapolations, which give 

m^2 l = 5.613(66) (MeV), (V.68) 
mf Gl {K) = 147.1(17) (MeV), (V.69) 
mf G1 ((j)) = 187.1(41) (MeV). (V.70) 

To compare the previous results, these values are converted to the MS scheme by 



where four loop expression is used for renormalization group functions f3 and r in the MS scheme 30, 



38]. The results in Fig. 1211 show good scalings, and the constant continuum extrapolations give 

m^f (NDR, 2GeV) = 4.026(48) (MeV), (V.72) 

mf^(ET) (NDR,2GeV) = 105.6(12) (MeV), (V.73) 

mf^(0)(NDR,2GeV) = 134.2(30) (MeV). (V.74) 

Contrary to the case of Bk, perturbatively renormalized quark masses of the previous CP-PACS 
result (filled squares) are underestimated, as seen in the figure. This clearly shows that the necessity 
of the non-perturbative renormalization for precision calculations in lattice QCD. We think that 
the large effects of the renormalizations are mainly canceled in the ratio of the Bk definition, and 
therefore, such cancellations can not be expected for general operators. 



VI. CONCLUSION AND DISCUSSIONS 

In this paper we have performed the non-perturbative renormalization of Bk and quark masses 
in the quenched domain-wall QCD using the Schrodinger functional method. Combined with the 
non-perturbative running obtained by Alpha collaboration, we have obtained the renormalization 
factors, which convert the lattice bare Bk and quark masses previously obtained by the CP-PACS 
collaboration to the RGI values. We obtain 

B K = 0.782(5)(7), (VI.75) 
mRGi = 5 . 613 ( 66 ) (MeV), (VI.76) 

mf GI {K) = 147.1(17) (MeV), (VI.77) 
rhf GI {(/)) = 187.1(41) (MeV). (VI.78) 
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in the continuum limit. These values correspond to renormalized values in MS scheme with the 
naive dimensional regularization are given as 

£]P(NDR, 2GeV) = 0.565(4)(5), (VI.79) 

m^(NDR, 2GeV) = 4.026(48) (MeV), (VI.80) 

mf^iK) (NDR,2GeV) = 105.6(12) (MeV), (VI.81) 

mp(</>)(NDR,2GeV) = 134.2(30) (MeV). (VI.82) 

With the non-perturbative renormalization in the DWQCD and data at a -1 ~ 2,3,4 GeV, we 
can extract Bk within 2% errors, except the quenching errors. The error in Bk is directly reflected 
to that in the CKM triangle constraint from ex- If we adopt our Bk for an input, the error of the 
constraint is improved as is shown in Fig. [22l The solid lines are central value and one standard 
deviation of the constraint from ex with our Bk- The dashed lines are results with Bk adopted 
by the CKM fitter group [42] • For other inputs we used those given by the CKM fitter. ' 

Although the perturbative renormalization can also achieve the same level of accuracy for Bk, 
it is clearly shown that the nonperturbative renormalization is indeed necessary for the precise 
determination of quark masses in lattice QCD. 
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APPENDIX A: CHIRAL WARD-TAKAHASHI IDENTITY FOR SF FORMALISM 
WITH DWF 



In this appendix we derive the Ward-Takahashi identity (|IV.58p for the Schrodinger functional 
formalism with domain-wall fermion. An explicit form of the chirally rotated correlation function 
used for evaluation of Zyv+AA is presented. We use the same notation of Ref. [13] for the SF 
formalism in this appendix. 

In this paper we adopt following massless DWF action with an orbifolding projection 



- l± loTo PQR 
ii ± = ^ • 



(A.l) 
(A.2) 



Here P is a parity transformation in fifth direction P Sj tifj(x, xo, t) = ip(x, xq, N5 — s + 1), and R is 
a time reflection operator acting on the temporal direction R XOtyo ip(x, yo, s) = ip(x, —xq,s). Q is 
the vector charge matrix for the chiral transformation 29[] 



Q 



s.t 



(\ \ 

1 

1 

0-10 

-1 

\0 -1/ 



(foriVs = 6). 



(A.3) 



For this action the physical quark propagator is given by 



^Lk(^>y) = 2(G quark n + ) x , y , n± 



l±7o« 



(A.4) 



This propagator is shown to agree with ()III.37[) numerically and we employ the latter in our 
numerical simulations. 

We consider the chiral rotation of the first flavour 



tfi(a;,a)-^0r,a) = {e~ iaQ ) st ih(x, t), ^(x, s) -> s) = $ x {x, t) {e iaQ ) ts , 
under which the physical quark and boundary quark fields are rotated as 



Qi(x), qi(x)=q 1 (x)e 



C[(x) = e-^CHx), Ux) = C[(x)e-^. 
For a = tt/2 the action for the first flavour is transformed as 

cSF _^oSF , v 
D dwf — * D dwf ' 1 i 

Y = a 4 ^#T_X^, 

aX = [PL5 s N^5 t N^ +l + PR8 s jf£ +1 8 t ^) S X)V . 



(A.5) 

(A.6) 
(A.7) 
(A.8) 

(A.9) 
(A.10) 

(A.ll) 
(A.12) 
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We notice that the orbifolding projection in the rotated action has an opposite sign and the 
rotated quark fields satisfy the opposite Dirichlet boundary condition to (|III.39|) 

P_§1 (x) \ xo=0 = 0, P + qi(x)\ X0=NT = 0, (A.13) 
%x( x ) p +\xo=o = °> ?i(x)P-\ X0 =N T = 0. (A.14) 

Now we derive the Ward-Takahashi identity for the four point function (|II.18|) and two point 
functions (|II.23|) as 



^aTbTo = ^2i[r A ]O i5 [T B ]0$ A+AV (x)0' 53 [T c i)s 

= ^(d2i[T A il 5 }0 A5 [T B }id$ v+AA (x)0' 53 [Tc})s +Y , (A.15) 

fi = -^6(Oi2[75]0 2 i[75])s = -^(^i 2 [l]^2i[l])s + y, (A.16) 

1 3 1 3 ~ 

-^J2(°'M°^})s = - g^ ^(0[ 2 [iw k }0 21 [ lk i l5 ]) §+Y , (A.17) 

k=X k=l 



k 



where 

Ovv+aa = (qilnQi) (toWi) + (917^7592) (q 3 7^7594) 

± { (qiIvQ*) {qz!^2) + (qii^im) (937^75^2)} • (A.18) 

Operators with tilde consist of the chirally rotated field for the first flavour, which satisfies the 
opposite SF Dirichlet boundary condition. Subscript S and S + Y mean the action under which 
the VEV is taken. 

The chiral symmetry breaking effect comes from a contribution of Y on the right hand side of the 
WT identity, which generates an operator mixing with the parity even operator. The contribution 
should be suppressed exponentially in the physical quark propagator and this is the case at tree 
level. We evaluate the effect by comparing two correlation functions directly 

(0 21 [T A ]0 45 [T B ]0$ A+AV (x)0' 53 [T c ])s ~ {0 21 [T A i lb \0^[T B }id^ v+AA (x)0', 3 [T c ])s. (A.19) 
For this purpose we define the renormalization factor for the parity even four fermi operator as 

Z$V + AA;s(90, an) = C^To^ ' (A.20) 
h s [xq = L/2;g ) 

ht(xo) = ^%j} X °\ (A.21) 
fx 

*fW (*o) = ^{0 2l [T A i^\0^[T B ]id^ v+AA {x)0^[T c \)^ (A.22) 
1 

2T 6 



/i = -^76(^ 2 [l]^2i[% (A.23) 



and compare it with that for the parity odd operator. 



APPENDIX B: NUMERICAL ANALYSIS OF B K FOR = 3.2 



Since the numerical data of CP-PACS collaboration at /3 = 3.2 is new and not published, we 
^ive a short summary of its numerical analysis in this appendix. 
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1. Run parameters and measurements 

We carry out run at j3 = 3.2, corresponding to a lattice spacing a~ l = 4.19(6) GeV determined 
from the the p meson mass m p = 770 MeV. We use the lattice size N% xNtXN& = 48 3 x 80 x 16. This 
lattice has a reasonably large spatial size of aN a ~ 2.3 fm and fifth dimensional length N$ = 16, 
which has been confirmed to be enough for Bk 0]- In this numerical simulation the domain wall 
height is taken to be M = 1.8. 

We take degenerate quarks in our calculations. The value of bare quark mass is chosen to be 
m f a = 0.009,0.018,0.027,0.036, which covers the range that m PS /m v ps 0.5334 - 0.8224. 

Quenched gauge configurations are generated on four-dimensional lattices. A sweep of gauge 
update contains one pseudo-heatbath and four overrelaxation steps. After a thermalization of 2000 
sweeps hadron propagators and 3-point functions necessary to evaluate Bk are calculated at every 
200th sweep. The gauge configuration on each fifth dimensional coordinate s is identical and is 
fixed to the Coulomb gauge. 

The domain- wall quark propagator needed to extract the Bk is calculated by the conjugate 
gradient algorithm with an even-odd pre-conditioning. Two quark propagators are evaluated for 
each configuration corresponding to the wall sources placed at either t = 4 or 77 in the time 
direction with the Dirichlet boundary condition, while the periodic boundary condition is imposed 
in the spatial directions. The two quark propagators are combined to form the kaon Green function 
with an insertion of the four-quark operator at time slices 1 < t < Nt in a standard manner. We 
employ the same quark propagators to evaluate pseudoscalar and vector meson propagators, and 
extract their masses. 



2. Pseudo scalar and vector meson masses 



We extract pseudoscalar and vector meson masses mps and my at each mj by a single expo- 
nential fit with meson propagators. Representative plots of effective masses are shown in Figs. [231 
and 1241 Fitting ranges chosen from inspection of such plots are 24 < t < 55 and 14 < t < 30 for 
pseudoscalar and vector meson masses, respectively. 

For the chiral extrapolation we fit mp s and my linearly in mja as illustrated in Fig. [25j Since 
pseudoscalar meson mass does not vanish at m/ = 0, we employ a fit of the form 

m 2 PS a 2 = Aps(rrifa + m res a), (B.l) 
mya = Ay + Byrrifa (B.2) 

and determine the parameters Aps,m res a for the pseudoscalar meson, and Ay, By for the vector 
meson. The physical bare masses, for the generated u and d quarks and m s j for the s quark 
are determined by equations that 



Aps{mfa + m res a) 
Ay + Bymfa 



Aps ( (m s jr(K)a + m u j d a)/2 + m ves a 



Ay + B v m u f d a 



f 

Ay + B v m s f (cl>)a 
Ay + B v m u f d a 



m v _ 


0.135 


m p 


0.77 ' 


m K 


0.495 


m p 


0.77 ' 


m<f> __ 


1.0194 


m p 


0.77 



(B.3) 

(B.4) 
(B.5) 
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For the s quark, we extract two values of the s quark mass, m S j(K) from the kaon mass input or 
m s jr((p) from the phi meson mass input. We then fix the lattice spacing a by setting the vector 
meson mass at the physical quark mass point rrij to the experimental value ra p = 770 MeV. 
Numerical values of lattice spacing and quark masses are listed in table [I] 

3. Extraction of B parameters 

In the course of our simulation we measure the kaon Bk P-2p . The s and d quark fields defining 
Bk are the physical fields given by (|III.35j) . and the four-quark and bilinear operators are taken 
to be local in the 4-dimensional space-time. 

In Fig. [26] we show typical data for the ratio of kaon Green functions for Bk as a function of 
the temporal site t of the weak operator. The values of these quantities at each raj are extracted 
by fitting the plateau with a constant. The fitting range, determined by the inspection of plots for 
the ratio and those for the effective pseudo scalar meson mass, is 24 < t < 55. 

The bare value of Bk is interpolated as a function of raja using a formula suggested by chiral 
perturbation theory, 

Bk = B (1 — 3c ra^a log(mja) + b raja) . (B.6) 

This interpolation is illustrated in Fig. [271 The physical value of Bk is obtained at the point 
rrij = (m S j{K)+m U j d ) /2 estimated from the experimental value of tuk /rap (a solid circle in Fig. [27]). 
The result is given in table E 
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2.6 2.9 3.2 


B (o) 


0.5908(57) 0.5655(69) 0.5478(71) 


a- 1 (GeV) 


1.807(37) 2.807(55) 4.186(65) 


m ud (MeV) 
m ud + m rcs (MeV) 


0.40(34) 1.85(50) 4.65(52) 
3.306(68) 3.226(72) 3.216(52) 


m s (K) (MeV) 
m s (K)+m rcs (MeV) 
m s ((f>) 


83.8(19) 83.2(20) 85.8(14) 
86.7(18) 84.6(19) 84.3(14) 
114.1(93) (MeV) 104.6(50) 107.8(27) 



TABLE I: Numerical values of bare B K in DWQCD given by CP-PACS collaboration [6j. Values of lattice 
spacing and bare quark masses are also listed. Here m u d is the u, d quark mass from the ir input, m s (K / <j)) 
is the strange quark mass from the K (0) input, and Tribes IS 0< residual quark mass where the pion mass 
vanishes. Data at = 3.2 is new and not published in Ref. p]. 



N L 


6 


8 


10 


12 


14 


16 


18 


20 




2.4446 


2.6339 


2.7873 


2.9175 


3.0313 


3.1331 


3.2254 


3.3103 


num. of conf. for scheme 1^5 
num. of conf. for scheme 6^9 


5165 
1165 


3632 
1032 


2000 
1000 


2188 
670 


1000 

284 


868 
312 


778 
200 


104 
104 



TABLE II: Values of (3 which satisfies clNt = 2L max and number of configurations for each lattice size. Data 
at Nl — 20 is used only for the step scaling function. 
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Z VA 

Z 



VA 



z 



VA 



^VA 
Z VA 
Z VA 



VA 



^VA 

7+ 
^VA 



J VA 



VA 



VA 



VA 



J VA 



J VA 



VA 



VA 



VA 



Zb k 
Zb k 
z b k 
z b k 
Zb k 
Zb k 
Zb k 
Zb k 
Zb k 



Z P 
Z m 



6 

2.4446 



2.6339 



_L0 

2.7873 



_L2 

2.9175 



14 

3.0313 



_L6 

3.1331 



_L8 

3.2254 



0.79580(23) 0.85837(18) 0.87048(15) 0.90395(12) 0.92286(30) 0.94463(13) 0.96184(12) 



0.69767(89) 

0.7196(13) 

0.7398(11) 

0.6850(11) 

0.6819(11) 

0.6382(18) 

0.6817(16) 

0.6319(17) 

0.6292(18) 



0.8778(14 
0.9089(21 
0.9439(18 
0.8528(17 
0.8489(18 
0.7812(25 
0.8541(21 
0.7702(23 
0.7668(25 



0.61284(68) 

0.6677(14) 

0.69325(88) 

0.6082(11) 

0.6074(11) 

0.5924(21) 

0.6391(13) 

0.5614(19) 

0.5608(19) 



0.7101(11 
0.7791(22 
0.8290(15 
0.6931(17 
0.6910(17 
0.6679(28 
0.7506(20 
0.6245(25 
0.6233(25 



1.1017(12) 
1.1363(18) 
1.1682(15) 
1.0816(15) 
1.0768(15) 
1.0084(23) 
1.0772(20) 
0.9984(22) 
0.9942(23) 



1.1914(19 
1.2336(28 
1.2811(24 
1.1574(23 
1.1521(23 
1.0605(32 
1.1594(27 
1.0456(30 
1.0410(32 



0.9189(18 
0.9679(28 
0.9956(23 
0.9030(23 
0.8985(23 
0.8224(26 
0.8930(23 
0.8103(24 
0.8061(26 



0.7024(13 
0.7937(28 
0.8301(18 
0.6974(21 
0.6968(21 
0.6726(27 
0.7433(19 
0.6244(24 
0.6235(24 



1.2127(24 
1.2774(37 
1.3139(30 
1.1918(29 
1.1857(30 
1.0854(34 
1.1787(30 
1.0695(31 
1.0640(33 



1.0179(21 
1.0758(31 
1.1088(26 
0.9982(24 
0.9933(25 
0.8980(35 
0.9822(31 
0.8838(32 
0.8785(34 



0.7472(15 
0.8516(28 
0.8945(20 
0.7413(21 
0.7400(20 
0.7135(36 
0.7958(27 
0.6572(32 
0.6570(32 



1.2457(25 
1.3166(38 
1.3570(32 
1.2216(30 
1.2157(30 
1.0990(42 
1.2021(37 
1.0817(39 
1.0752(41 



0.65512(64) 
1.5264(15) 



0.66259(88) 
1.5092(20) 



0.64990(99) 
1.5387(23) 



0.6560(10 
1.5244(24 



1.0785(34 
1.1521(53 
1.1804(44 
1.0648(40 
1.0595(42 
0.9592(53 
1.0450(49 
0.9422(48 
0.9395(51 



0.7587(25 
0.8820(48 
0.9152(33 
0.7630(35 
0.7627(34 
0.7295(55 
0.8089(44 
0.6698(48 
0.6706(48 



1.2663(40 
1.3527(63 
1.3860(52 
1.2502(48 
1.2441(50 
1.1285(62 
1.2294(57 
1.1086(56 
1.1054(60 



0.6580(16 
1.5198(37 



1.1388(37 
1.2132(59 
1.2484(48 
1.1196(45 
1.1138(47 
1.0093(56 
1.1058(49 
0.9913(50 
0.9859(55 



0.7886(25 
0.9074(53 
0.9532(35 
0.7837(38 
0.7823(38 
0.7537(56 
0.8418(44 
0.6921(50 
0.6909(50 



1.2763(42 
1.3596(66 
1.3991(55 
1.2547(50 
1.2482(53 
1.1313(62 
1.2394(55 
1.1111(56 
1.1050(61 



0.6606(16 
1.5137(37 



1.1821(44) 

1.264(18) 

1.2994(58) 

1.1645(54) 

1.157(16) 

1.039(31) 

1.1409(66) 

1.0202(65) 

1.015(28) 



0.8061(30) 

0.934(13) 

0.9783(41) 

0.8019(42) 

0.802(10) 

0.769(23) 

0.8618(55) 

0.7018(55) 

0.703(19) 



1.2778(47) 

1.366(19) 

1.4045(63) 

1.2587(58) 

1.251(17) 

1.123(33) 

1.2329(72) 

1.1025(71) 

1.097(30) 



0.6627(18) 
1.5091(42) 



TABLE III: Numerical values of Z v {ga), Zy A+Av . s (g ,afj, m in) and Z BK - s {go, »Mmin) at 2L max . Values of 
Zp(go 7 a/imin) and Z m (go, a/j, m i n ) are also listed for the latter use. 
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N L 


4 


6 


8 


10 


20 


P 


2.6339 


2.9175 


3.1331 


3.3103 


3.3103 


Z v 


0.9905(21) 


0.87148(98) 


0.97930(77) 


0.97587(57) 


0.97898(26) 


Z VA;1 


1.0784(69) 


0.7937(48) 


1.0967(59) 


1.0853(70) 


1.258(11) 


7+ 


1.0302(74) 


0.8061(63) 


1.0945(84) 


1.1054(91) 


1.349(16) 


7+ 


1.1000(72) 


0.8215(56) 


1.1446(67) 


1.1384(83) 


1.376(14) 


7+ 
VAA 


1.0058(70) 


0.7807(56) 


1.0480(74) 


1.0555(76) 


1.246(12) 


7+ 


1.0076(71) 


0.7811(57) 


1.0498(77) 


1.0555(76) 


1.240(12) 


7+ 


1.0018(70) 


0.7703(54) 


1.0346(73) 


1.0366(73) 


1.1181(91) 


7+ 

^VA-,7 


1.0797(69) 


0.7969(49) 


1.1024(59) 


1.0907(71) 


1.2145(89) 


7+ 

^VA-,8 


0.9871(67) 


0.7574(51) 


1.0094(68) 


1.0113(66) 


1.0995(85) 


7+ 

^VA-,9 


0.9890(68) 


0.7578(51) 


1.0112(72) 


1.0112(66) 


1.0942(88) 


Z VA;1 


1.0609(63) 


0.7389(30) 


0.9721(49) 


0.9349(57) 


0.8252(67) 


Z VA;2 


0.9831(82) 


0.7746(69) 


0.9775(93) 


0.9847(99) 


0.971(15) 


Z VA : 3 


1.1135(69) 


0.7933(42) 


1.0644(60) 


1.0334(74) 


0.995(10) 


Z VA;4 


0.9423(74) 


0.7288(58) 


0.9043(79) 


0.9033(79) 


0.842(11) 


Z VA;5 


0.9435(74) 


0.7329(57) 


0.9088(79) 


0.9070(75) 


0.841(11) 


Z VA;6 


0.9559(77) 


0.7401(61) 


0.9239(83) 


0.9234(81) 


0.805(11) 


Z VA;7 


1.0929(66) 


0.7696(37) 


1.0252(54) 


0.9901(65) 


0.8777(72) 


Z VA;8 


0.9248(71) 


0.7070(53) 


0.8710(73) 


0.8654(70) 


0.7430(93) 


Z VA;9 


0.9260(71) 


0.7110(54) 


0.8753(74) 


0.8690(67) 


0.7425(93) 


Zb k ;1 


1.0992(51) 


1.0450(53) 


1.1436(61) 


1.1397(73) 


1.313(11) 


Zb k ;2 


1.0502(59) 


1.0615(74) 


1.1413(88) 


1.1608(96) 


1.407(17) 


Z Bk;3 


1.1213(55) 


1.0816(63) 


1.1935(71) 


1.1954(87) 


1.436(15) 


Zbk'A 


1.0252(54) 


1.0280(64) 


1.0928(76) 


1.1084(80) 


1.300(13) 


Z B k ;5 


1.0271(55) 


1.0285(65) 


1.0947(80) 


1.1084(81) 


1.294(13) 


Zb k ;G 


1.0211(54) 


1.0142(62) 


1.0788(75) 


1.0885(76) 


1.1667(95) 


Z B K ;7 


1.1005(51) 


1.0493(54) 


1.1495(62) 


1.1453(74) 


1.2672(93) 


Z B K :8 


1.0062(51) 


0.9972(57) 


1.0525(69) 


1.0619(69) 


1.1472(89) 


Z B K ;9 


1.0081(52) 


0.9977(58) 


1.0544(73) 


1.0619(70) 


1.1417(93) 


z P 


0.8718(36) 


0.7680(33) 


0.8160(36) 


0.8051(37) 


0.6712(53) 




1.1471(48) 


1.3020(55) 


1.2254(54) 


1.2421(56) 


1.490(12) 



TABLE IV: Numerical values of Z v (g ), Zy A+Av . s {g ,2a^ m - m ), Z BK;s (ga, 2a^ min ), Z P (g , 2a/j, mhl ) and 
Z m (go, 2a/z m in) at L max for SSF. A box size of Nl — 20 lattice is 2L max . 
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4 


6 


8 


10 


8 


12 


16 


20 


P 


2.6339 


2.9175 


3.1331 


3.3103 


2.6339 


2.9175 


3.1331 


3.3103 


7+ 


0.7664(38) 


0.7669(49) 


0.7928(54) 


0.8157(50) 


0.8778(14) 


1.0179(21) 


1.1388(37) 


1.258(11) 


7+ 


0.7602(44) 


0.7661(60) 


0.7977(71) 


0.8210(71) 


0.9089(21) 


1.0758(31) 


1.2132(59) 


1.349(16) 


7+ 


0.7829(39) 


0.7937(55) 


0.8260(62) 


0.8629(96) 


0.9439(18) 


1.1088(26) 


1.2484(48) 


1.376(14) 


7+ 

^VA-A 


0.7440(42) 


0.7397(52) 


0.7655(60) 


0.801(14) 


0.8528(17) 


0.9982(24) 


1.1196(45) 


1.246(12) 


Z VA;5 


0.7445(42) 


0.7407(52) 


0.7670(63) 


0.802(14) 


0.8489(18) 


0.9933(25) 


1.1138(47) 


1.240(12) 


7+ 


0.7379(41) 


0.7321(51) 


0.7542(58) 


0.7748(61) 


0.7812(25) 


0.8980(35) 


1.0093(56) 


1.1181(91) 


7+ 


0.7675(38) 


0.7700(50) 


0.7957(53) 


0.8302(89) 


0.8541(21) 


0.9822(31) 


1.1058(49) 


1.2145(89) 


7+ 

^VA-.S 


0.7294(40) 


0.7176(47) 


0.7374(52) 


0.771(13) 


0.7702(23) 


0.8838(32) 


0.9913(50) 


1.0995(85) 


7+ 
VA;9 


0.7298(41) 


0.7186(48) 


0.7389(56) 


0.771(14) 


0.7668(25) 


0.8785(34) 


0.9859(55) 


1.0942(88) 


ZvA:l 


0.8065(28) 


0.7611(37) 


0.7532(41) 


0.7437(39) 


0.7101(11) 


0.7472(15) 


0.7886(25) 


0.8252(67) 


ZvA:2 


0.8076(47) 


0.7693(60) 


0.7688(68) 


0.7573(81) 


0.7791(22) 


0.8516(28) 


0.9074(53) 


0.971(15) 


Z VA:3 


0.8465(34) 


0.8193(48) 


0.8221(52) 


0.8241(92) 


0.8290(15) 


0.8945(20) 


0.9532(35) 


0.995(10) 


Z VA;4 


0.7757(43) 


0.7211(49) 


0.7114(56) 


0.711(13) 


0.6931(17) 


0.7413(21) 


0.7837(38) 


0.842(11) 


Z VA:5 


0.7770(42) 


0.7247(47) 


0.7165(53) 


0.717(13) 


0.6910(17) 


0.7400(20) 


0.7823(38) 


0.841(11) 


Z VA:6 


0.7839(44) 


0.7352(51) 


0.7269(56) 


0.7147(70) 


0.6679(28) 


0.7135(36) 


0.7537(56) 


0.805(11) 


Z VA:7 


0.8299(32) 


0.7948(44) 


0.7920(45) 


0.7929(85) 


0.7506(20) 


0.7958(27) 


0.8418(44) 


0.8777(72) 


Z VA:& 


0.7605(41) 


0.6996(45) 


0.6854(50) 


0.684(12) 


0.6245(25) 


0.6572(32) 


0.6921(50) 


0.7430(93) 


Z VA\9 


0.7617(40) 


0.7031(43) 


0.6902(48) 


0.690(12) 


0.6233(25) 


0.6570(32) 


0.6909(50) 


0.7425(93) 


Zb k ;1 


1.0096(33) 


1.0715(63) 


1.1041(69) 


1.1221(66) 


1.1914(19) 


1.2457(25) 


1.2763(42) 


1.313(11) 


Zb k \2 


1.0014(41) 


1.0704(79) 


1.1110(94) 


1.1295(96) 


1.2336(28) 


1.3166(38) 


1.3596(66) 


1.407(17) 


Zb k ;3 


1.0314(36) 


1.1089(72) 


1.1504(81) 


1.187(13) 


1.2811(24) 


1.3570(32) 


1.3991(55) 


1.436(15) 


Zb k -4 


0.9801(38) 


1.0334(67) 


1.0661(78) 


1.102(19) 


1.1574(23) 


1.2216(30) 


1.2547(50) 


1.300(13) 


Z B K ;5 


0.9807(39) 


1.0348(68) 


1.0683(83) 


1.103(19) 


1.1521(23) 


1.2157(30) 


1.2482(53) 


1.294(13) 


Z B K ;6 


0.9721(37) 


1.0228(65) 


1.0504(76) 


1.0659(82) 


1.0605(32) 


1.0990(42) 


1.1313(62) 


1.1667(95) 


ZbkJ 


1.0112(33) 


1.0758(63) 


1.1082(68) 


1.142(12) 


1.1594(27) 


1.2021(37) 


1.2394(55) 


1.2672(93) 


Z B K :» 


0.9608(36) 


1.0026(59) 


1.0270(68) 


1.060(18) 


1.0456(30) 


1.0817(39) 


1.1111(56) 


1.1472(89) 


Zb k :9 


0.9615(37) 


1.0040(60) 


1.0291(73) 


1.061(18) 


1.0410(32) 


1.0752(41) 


1.1050(61) 


1.1417(93) 


z P 


0.8277(22) 


0.7670(30) 


0.7394(35) 


0.7240(31) 


0.66259(88) 


0.6560(10) 


0.6606(16) 


0.6712(53) 


Z m 


1.2082(32) 


1.3038(51) 


1.3524(65) 


1.3812(60) 


1.5092(20) 


1.5244(24) 


1.5137(37) 


1.490(12) 



TABLE V: Numerical values of Zy A+Av . s (g ,a/j,), Z BK ; 8 (go,afj), Z P (g ,a/j,) and Z m {g Q ,a^) at L max and 
2L max at M — 1.4 for SSF. 
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7+ 
^VA 

Z VA 
7+ 

7+ 
^VA 

Z VA 

7+ 
^VA 

Z VA 

7+ 
^VA 

7+ 
^VA 



J VA 



J VA 



J VA 



J VA 



J VA 



J VA 



J VA 



J VA 



z 



VA 



Zb k 
Zb k 
z b k 
z b k 
z b k 
Zb k 
Zb k 
Zb k 
Zb k 



Z P 
Z m 



6 

2.4446 



2.6339 



10 

2.7873 



_L2 

2.9175 



U 

3.0313 



_L6 

3.1331 



_L8 

3.2254 



0.85609(24) 0.88162(19) 0.88822(16) 0.91166(12) 0.92852(30) 0.94719(13) 0.96359(12) 



0.8285(11 
0.8546(15 
0.8785(13 
0.8134(13 
0.8098(13 
0.7582(21 
0.8099(18 
0.7507(20 
0.7475(21 



0.72775(80) 

0.7929(17 

0.8232(10 

0.7223(13 

0.7213(13 

0.7038(24 

0.7593(16 

0.6669(22 

0.6663(22 



1.1304(12 
1.1660(18 
1.1987(15 
1.1098(15 
1.1049(16 
1.0351(24 
1.1056(20 
1.0248(22 
1.0204(24 



0.9344(15 
0.9675(23 
1.0047(19 
0.9077(18 
0.9036(19 
0.8317(26 
0.9093(23 
0.8200(25 
0.8164(26 



0.7558(12 
0.8293(23 
0.8825(16 
0.7377(18 
0.7355(18 
0.7111(30 
0.7991(21 
0.6649(26 
0.6636(27 



1.2021(19 
1.2447(28 
1.2927(24 
1.1679(23 
1.1626(23 
1.0703(32 
1.1701(27 
1.0553(30 
1.0506(32 



0.70476(69) 
1.4189(14) 



0.68054(9 
1.4694(20 



TABLE VI: Numerical values of Z v 



0.9645(19 
1.0160(30 
1.0450(24 
0.9479(24 
0.9431(25 
0.8634(28 
0.9375(24 
0.8507(26 
0.8463(27 



0.7373(14 
0.8331(29 
0.8713(19 
0.7320(22 
0.7314(22 
0.7061(29 
0.7803(20 
0.6555(25 
0.6546(26 



1.2225(24 
1.2878(37 
1.3246(31 
1.2015(30 
1.1954(31 
1.0945(34 
1.1885(30 
1.0784(31 
1.0729(33 



0.6631(10 
1.5080(23 



1.0392(21 
1.0983(32 
1.1320(27 
1.0191(25 
1.0141(26 
0.9169(35 
1.0029(31 
0.9024(33 
0.8970(35 



0.7628(16 
0.8694(29 
0.9132(20 
0.7568(22 
0.7555(21 
0.7285(36 
0.8126(28 
0.6710(33 
0.6708(33 



1.2504(25 
1.3215(38 
1.3621(32 
1.2262(30 
1.2202(31 
1.1032(42 
1.2067(37 
1.0859(39 
1.0794(41 



0.6616(10 
1.5115(24 



1.0951(34 
1.1698(53 
1.1985(44 
1.0811(41 
1.0758(42 
0.9741(53 
1.0612(50 
0.9568(48 
0.9541(52 



0.7704(25 
0.8956(48 
0.9293(33 
0.7748(35 
0.7745(34 
0.7408(56 
0.8214(45 
0.6802(49 
0.6811(49 



1.2702(40 
1.3568(63 
1.3902(52 
1.2540(48 
1.2479(50 
1.1321(62 
1.2334(58 
1.1121(56 
1.1089(60 



0.6620(16 
1.5106(37 



1.1470 
1.2218 
1.2573 
1.1276 
1.1218 
1.0166 
1.1138 
0.9985 
0.9930 



0.7942 
0.9139 
0.9600 
0.7893 
0.7879 
0.7592 
0.8479 
0.6971 
0.6959 



1.2784 
1.3619 
1.4014 
1.2568 
1.2503 
1.1333 
1.2416 
1.1131 
1.1070 



0.6624 
1.5096 



1.1880(44) 

1.270(18) 

1.3058(58) 

1.1703(54) 

1.163(16) 

1.044(31) 

1.1467(66) 

1.0254(66) 

1.020(28) 



0.8101(30) 

0.938(13) 

0.9832(41) 

0.8058(43) 

0.806(10) 

0.773(23) 

0.8662(55) 

0.7054(56) 

0.706(19) 



1.2795(47) 

1.368(19) 

1.4064(63) 

1.2604(59) 

1.253(17) 

1.124(33) 

1.2347(72) 

1.1041(71) 

1.099(31) 



0.6639(18) 
1.5063(42) 



go), Zy A+Av . s {go, a/Xmin) and Z Bk ; S {go, a/i m in) at 2L max with the tree 



level improved renormalization condition (jIV.43|) . Values of Zp(go, a/i m i n ) and Z m (go, o/i m i n ) are also listed. 
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N L 


4 


6 


8 


10 


20 


P 


2.6339 


2.9175 


3.1331 


3.3103 


3.3103 


Z v 


1.1071(24) 


0.9756(11) 


1.01320(80) 


1.00782(59) 


0.97970(26) 


7+ 


1.3728(88) 


1.0330(63) 


1.1844(64) 


1.1716(76) 


1.261(11) 


7+ 

^VA-,2 


1.3115(94) 


1.0493(83) 


1.1820(91) 


1.1933(98) 


1.352(16) 


7+ 

^VA-3 


1.4004(92) 


1.0693(72) 


1.2361(72) 


1.2289(90) 


1.380(14) 


7+ 


1.2803(89) 


1.0162(73) 


1.1318(80) 


1.1395(82) 


1.249(12) 


7+ 


1.2827(90) 


1.0167(74) 


1.1338(84) 


1.1394(82) 


1.243(12) 


7+ 

^VA-fi 


1.2753(89) 


1.0026(71) 


1.1173(79) 


1.1190(79) 


1.1209(91) 


7+ 


1.3744(88) 


1.0373(64) 


1.1906(64) 


1.1774(77) 


1.2175(89) 


7+ 


1.2566(85) 


0.9858(66) 


1.0901(74) 


1.0917(71) 


1.1022(86) 


7+ 

^VA-9 


1.2590(87) 


0.9863(67) 


1.0920(77) 


1.0916(72) 


1.0968(89) 


Z VA;1 


1.3505(80) 


0.9618(40) 


1.0498(53) 


1.0092(62) 


0.8273(68) 


Z VA;2 


1.251(10) 


1.0083(90) 


1.056(10) 


1.063(11) 


0.973(15) 


Z VA-3 


1.4175(88) 


1.0325(55) 


1.1496(65) 


1.1156(79) 


0.997(10) 


Z VA;4 


1.1995(94) 


0.9486(75) 


0.9766(85) 


0.9751(85) 


0.844(11) 


Z VA-5 


1.2010(94) 


0.9539(75) 


0.9815(85) 


0.9791(81) 


0.843(11) 


Z VA-fi 


1.2169(98) 


0.9634(79) 


0.9978(89) 


0.9968(87) 


0.807(11) 


Z VA-7 


1.3913(84) 


1.0017(48) 


1.1072(58) 


1.0688(70) 


0.8799(73) 


Z VA;S 


1.1773(90) 


0.9203(69) 


0.9407(79) 


0.9342(76) 


0.7448(94) 


Z VA:9 


1.1788(90) 


0.9254(70) 


0.9453(80) 


0.9381(72) 


0.7443(93) 


Zb k ;1 


1.1200(52) 


1.0854(55) 


1.1537(62) 


1.1535(74) 


1.314(11) 


Zb k \2 


1.0700(60) 


1.1025(77) 


1.1514(89) 


1.1749(98) 


1.408(17) 


Z B k ;3 


1.1425(56) 


1.1234(65) 


1.2041(72) 


1.2099(88) 


1.437(15) 


Zbk'A 


1.0446(55) 


1.0677(67) 


1.1025(77) 


1.1218(81) 


1.301(13) 


Z B K ;5 


1.0465(56) 


1.0682(68) 


1.1044(81) 


1.1218(82) 


1.295(13) 


Z B K ;6 


1.0404(55) 


1.0534(64) 


1.0884(76) 


1.1017(77) 


1.1678(95) 


Zb k j 


1.1213(52) 


1.0898(56) 


1.1597(62) 


1.1592(75) 


1.2685(93) 


Zb k -?, 


1.0252(52) 


1.0358(60) 


1.0619(70) 


1.0748(70) 


1.1483(89) 


Z B K ;9 


1.0271(53) 


1.0363(60) 


1.0637(74) 


1.0748(70) 


1.1427(93) 


z P 


0.9745(41) 


0.8598(36) 


0.8443(37) 


0.8315(38) 


0.6717(53) 


Z m 


1.0262(43) 


1.1631(49) 


1.1844(52) 


1.2027(55) 


1.489(12) 



TABLE VII: Numerical values of Z v (go), z y A+Av . s (g ,2afi m in), Z BK ; S (go, 2a/imin), Zp(gg, 2a^ min ) and 
Z m (go, 2a/x m i n ) at L max for SSF with the tree level improved renormalization condition (|IV.43jl . A box size 
of Nl — 20 lattice is 2L max . 
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scheme 


°V A+AV:s 


a V A+AV\s 




1 


1.136(22) 


0.841(16) 


1.116(23) 


2 


1.221(33) 


0.988(32) 


1.197(34) 


3 


1.164(29) 


0.889(22) 


1.143(30) 


4 


1.182(30) 


0.944(28) 


1.158(31) 


5 


1.180(30) 


0.936(27) 


1.154(32) 


6 


1.120(24) 


0.905(26) 


1.099(25) 


7 


1.096(22) 


0.835(18) 


1.076(23) 


8 


1.113(24) 


0.887(25) 


1.091(25) 


9 


1.111(25) 


0.880(25) 


1.087(26) 



TABLE VIII: Value of the SSF in the continuum limit ^va+av-s an( ^ <7 b k -s at L max , obtained by a combined 
linear fit of data at M = 1.8 with the tree level improved renormalization condition and at M = 1.4. 



2.444602 2.633865 2.787275 2.917468 3.031335 3.133065 3.225406 



Z B K -1 
Z B K ;2 

Z B k ;4 

Zbk;5 
Z B K fi 

Zb k ;8 
Zb k ;9 



1.256(22) 
1.252(28) 
1.208(23) 
1.321(27) 
1.294(25) 
1.361(25) 
1.273(21) 
1.392(26) 
1.365(20) 



1.336(23) 
1.337(30) 
1.303(25) 
1.390(28) 
1.361(27) 
1.407(26) 
1.347(22) 
1.433(27) 
1.406(20) 



1.358(23) 
1.383(31) 
1.335(25) 
1.430(29) 
1.400(28) 
1.439(27) 
1.368(23) 
1.464(27) 
1.435(21) 



1.389(24) 
1.419(32) 
1.373(26) 
1.459(30) 
1.429(28) 
1.451(27) 
1.389(23) 
1.475(28) 
1.444(21) 



1.411(25) 
1.457(33) 
1.401(27) 
1.492(31) 
1.461(29) 
1.489(28) 
1.420(24) 
1.510(29) 
1.484(23) 



1.420(25) 
1.463(33) 
1.413(27) 
1.496(31) 
1.464(29) 
1.490(28) 
1.429(24) 
1.512(29) 
1.481(23) 



1.422(25) 
1.469(39) 
1.418(27) 
1.500(31) 
1.467(35) 
1.478(51) 
1.421(25) 
1.499(29) 
1.470(46) 



TABLE IX: Numerical values of renormalization factors ZB K - s {go) for the RGI Bk- 
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2.6 


2.9 


3.2 


continuum 


xVclof 




(5o) 


1.314(13) 


1.389(13) 


1.422(18) 






Zb k :2 


(.9o) 


1.321(17) 


1.419(18) 


1.470(26) 








(9o) 


1.279(14) 


1.373(14) 


1.417(20) 








(.90) 


1.378(16) 


1.459(16) 


1.502(22) 






Zg K; 5 


(50) 


1.350(15) 


1.429(16) 


1.470(23) 








(50) 


1.400(15) 


1.458(15) 


1.493(28) 








(50) 


1.327(13) 


1.395(13) 


1.425(18) 








(50) 


1.428(15) 


1.482(15) 


1.508(21) 






Zb k ;9 


(50) 


1.399(12) 


1.452(12) 


1.486(23) 








1 


0.777(11) 


0.786(12) 


0.779(14) 


0.7830(91) 


0.14 




2 


0.781(12) 


0.802(14) 


0.806(18) 


0.804(11) 


0.03 




3 


0.756(11) 


0.776(13) 


0.776(15) 


0.776(98) 


_ 




4 


0.814(12) 


0.825(14) 


0.823(16) 


0.824(11) 


0.009 


B K 


5 


0.797(12) 


0.808(13) 


0.805(17) 


0.807(10) 


0.019 


Bx 


6 


0.827(12) 


0.825(13) 


0.818(18) 


0.823(10) 


0.10 


Bx 


7 


0.784(11) 


0.789(12) 


0.781(14) 


0.7856(91) 


0.19 


Bx 


8 


0.844(12) 


0.838(13) 


0.826(16) 


0.833(10) 


0.34 


Bx 


9 


0.827(11) 


0.821(12) 


0.814(16) 


0.8185(96) 


0.12 



TABLE X: Renormalization factors ZB K -,s{go) and RGI Bk at three /3's with their continuum extrapolations 
for schemes s = 1, • • • ,9. Values of x 2 /dof are also listed. 



7MS 

7MS 

B. 

gMS 
Bk 
7MS 

Bk 

7MS 
7MS 



;l(ffO, 
; 2 (.90, 
53(50, 
54(50, 
;5(50, 
;6(50, 
;7(50, 
58 (50, 

; g(5o, 



2GeV) 
2GeV) 
2GeV) 
2GeV) 
2GeV) 
2GcV) 
2GeV) 
2GeV) 
2GeV) 



S||(NDR, 
S||(NDR, 
Sg(NDR, 
B^ S 5 (NDR, 
B|j(NDR, 
S|^(NDR, 
S|j(NDR, 

dMS 
n K:9 



(NDR. 



2GeV) 
2GeV) 
2GeV) 
2GeV) 
2GeV) 
2GeV) 
2GeV) 
2GeV) 
2GeV) 



2.6 



2.9 



3.2 



0.9494(96) 
0.954(12) 
0.924(10) 
0.995(12) 
0.975(11) 
1.011(11) 
0.9582(94) 
1.032(11) 
1.0108(90) 



1.003(10) 
1.025(13) 
0.991(11) 
1.054(12) 
1.032(12) 
1.053(11) 
1.0080(98) 
1.070(12) 
1.0490(91) 



1.027(13) 
1.062(19) 
1.023(15) 
1.085(16) 
1.062(17) 
1.078(20) 
1.029(13) 
1.089(16) 
1.073(17) 



0.5609(78) 
0.5638(92) 
0.5458(80) 
0.5880(90) 
0.5759(87) 
0.5974(88) 
0.5661(78) 
0.6095(89) 
0.5972(78) 



0.5675(89) 
0.580(10) 
0.5607(92) 
0.596(10) 
0.5837(97) 
0.5956(97) 
0.5700(89) 
0.6053(99) 
0.5932(89) 



0.563(10) 
0.582(13) 
0.561(11) 
0.594(12) 
0.582(12) 
0.591(14) 
0.564(10) 
0.597(12) 
0.588(12) 



continuum 



0.5655(66 
0.5807(79 
0.5608(71 
0.5952(77 
0.5830(75 
0.5941(80 
0.5673(66 
0.6019(76 
0.5914(71 



X 2 /dof 



0.11 
0.015 
0.0004 
0.016 
0.012 
0.073 
0.20 
0.28 
0.12 



TABLE XI: Renormalization factors Z^. s (g , 2GeV) and renormalized Bk in MS scheme with NDR at 
three (3's with their continuum extrapolations for schemes s = 1, • • • ,9. Values of % 2 /dof are also listed. 



28 



10 



12 



14 



16 



18 







2.444602 2.633865 2.787275 2.917468 3.031335 3.133065 3.225406 



7+ 

^VV+AA-,1 
7+ 

^VV+AA-2 

z+ 



z 



J VV+AA:3 
■+ 

VV+AA;4 
7+ 

7+ 

■^vv+AA-e 

7+ 

^VV+AA-J 
7+ 

^VV+AA-,8 
7+ 



0.7240(31 
0.7460(43 
0.7669(37 
0.7102(37 
0.7071(38 
0.6639(33 
0.7095(29 
0.6570(32 
0.6542(33 



Z 



J VV+AA;1 
7j VV+AA:2 
VV+AA:3 
7j VV+AA-A 
7 VV+AA;5 
7 VV+AA;6 
7j VA+AV-7 
VV+AA-fi 

7 vv+AA;9 



z 



0.6328(24 
0.6903(47 
0.7156(31 
0.6299(38 
0.6278(38 
0.6142(37 
0.6621(26 
0.5828(34 
0.5808(34 



0.8799(37 
0.9086(55 
0.9475(47 
0.8527(45 
0.8473(45 
0.7827(37 
0.8578(34 
0.7720(35 
0.7671(37 



0.7153(30 
0.7766(58 
0.8347(39 
0.6916(45 
0.6884(43 
0.6690(43 
0.7557(32 
0.6261(37 
0.6233(38 



0.9230(46 
0.9720(72 
1.0012(58 
0.9054(58 
0.9007(60 
0.8224(48 
0.8957(41 
0.8099(45 
0.8058(47 



0.7056(34 
0.7986(75 
0.8368(47 
0.7004(56 
0.6990(55 
0.6757(53 
0.7486(37 
0.6266(46 
0.6253(46 



1.0158(94) 

1.072(14) 

1.103(12) 

0.997(11) 

0.992(12) 

0.9061(96) 

0.9865(84) 

0.8912(88) 

0.8866(94) 



1.077(10) 

1.156(15) 

1.178(12) 

1.069(11) 

1.061(12) 

0.9581(89) 

1.0390(75) 

0.9431(80) 

0.9361(89) 



1.138(10) 

1.213(17) 

1.251(13) 

1.118(13) 

1.109(13) 

1.001(11) 

1.1008(91) 

0.9838(97) 

0.976(10) 



1.176(13) 

1.236(18) 

1.284(16) 

1.142(14) 

1.138(14) 

1.029(11) 

1.137(11) 

1.0103(99) 

1.008(11) 



0.7480(64) 

0.848(14) 

0.8842(91) 

0.741(11) 

0.743(11) 

0.716(11) 

0.7904(70) 

0.6626(94) 

0.6645(98) 



0.7614(73) 

0.912(14) 

0.9246(92) 

0.7914(96) 

0.7829(97) 

0.7559(92) 

0.8157(75) 

0.6982(76) 

0.6906(81) 



0.7928(86) 

0.909(15) 

0.964(11) 

0.787(11) 

0.781(11) 

0.750(10) 

0.8480(91) 

0.6928(94) 

0.6873(95) 



0.8108(94) 

0.907(14) 

0.975(13) 

0.7836(99) 

0.7856(99) 

0.7550(96) 

0.863(10) 

0.6935(85) 

0.6952(88) 



TABLE XII: Numerical values of renormalization factors Zy V+AA . s (g , a/i m ; n ) and Z vv+AA . s (g , a/j, m i n ) for 
the parity even operator with the chirally rotated scheme. A number of configurations is 100 on each lattice 



size. 





2.444602 2.633865 2.787275 2.917468 3.031335 3.133065 3.225406 


Z m 


1.642(17) 1.700(18) 1.745(18) 1.749(18) 1.748(19) 1.747(19) 1.743(19) 



TABLE XIII: Numerical values of the renormalization factor Z m (g ) for the RGI quark mass. 



/3 


2.6 2.9 3.2 


continuum 


X 2 /dof 


Z m 


1.694(10) 1.749(10) 1.742(13) 






m ud (MeV) 
m s (K) (MeV) 


0.68(58) 3.23(87) 8.11(90) 
141.9(33) 145.5(37) 149.4(27) 


12.2(15) 
154.8(52) 


2.0 
0.036 


m ud + m rcs (MeV) 
m s (K) +m ros (MeV) 
m s ((f>) (MeV) 


5.60(12) 5.64(13) 5.60(10) 
146.9(31) 147.9(34) 146.9(26) 
193(16) 183.0(88) 187.8(49) 


5.613(66) 
147.1(17) 
187.1(41) 


0.033 
0.033 
0.20 



TABLE XIV: The renormalization factor Z m (g ) and RGI light quark masses at three /3's with their con- 
tinuum linear extrapolations. Values of x 2 /dof are also listed. 
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p 


2.6 2.9 3.2 


continuum 


X 2 /dof 


Z™ s (MeV) 


1.2153(72) 1.2545(74) 1.2497(96) 






m^ s (MeV) 
mf s (K) (MeV) 


0.49(42) 2.31(63) 5.81(65) 
101.8(23) 104.4(26) 107.2(19) 


8.8(10) 
111.0(37) 


2.0 
0.036 


(m ud + m rcs ) MS (MeV) 
(m s (K) + m rcs ) MS (MeV) 
mf s (<f>) (MeV) 


4.018(86) 4.047(94) 4.020(72) 
105.3(23) 106.1(25) 105.4(19) 
139(11) 131.2(63) 134.7(35) 


4.026(48) 
105.6(12) 
134.2(30) 


0.033 
0.033 
0.20 



TABLE XV: The renormalization factor Z^ s (g ) and renormalized light quark masses in MS scheme at 
three /3's with their continuum linear extrapolations. Values of x 2 /dof are also listed. 
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L/a=6 



L/a=8 
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FIG. 1: x dependence of Zy(go) at various lattice sizes. 
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L7a=6 



L7a=8 



1.6 
1.2 
0.8 
0.4 
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FIG. 2: xq dependence of Zy j4+Ay . 1 (5 , a^min) for the scheme 1. 
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FIG. 3: x dependence of Zy A+Av . 2 (g 0} a/i m in) for the scheme 2. 
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FIG. 4: x dependence of Z yyl+Ay . 1 (5o, aMmin) for the scheme 1. 
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FIG. 5: xq dependence of Z VA+Av . 2 (g 0} a/i m in) for the scheme 2. 
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FIG. 7: The SSF T, VA+Av . 1 (u, a/L) as a function of the fifth dimensional length N5. Triangles (squares) are 
results at L/a — 4(6). 



36 



Tree level £v A+A vi 

(u,a/L) 




0.05 0.1 0.15 0.2 0.25 0.3 
a/L 



FIG. 8: Scaling behaviors of the tree level SSF. Open circles are results at M — 1.5 while open triangle at 
M = 0.9. The continuum value is represented by the dotted line. 
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FIG. 16: x dependence of Zp(g , l/2L max ) at various lattice sizes. 
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FIG. 17: The scaling behavior of the SSF T,p(u,a/L) with the ordinary renormalization condition. Open 
circles are our results with the domain- wall fermion at M — 1.8, while triangles show results by the Alpha 
collaboration with the improved Wilson fermion action. Corresponding continuum limits are represented by 
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FIG. 18: Scaling behaviors of the SSF £p(it, a/L) for the domain-wall fermion at M = 1.8 with the ordinary 
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FIG. 19: j3 dependences of the renormalization factor Z m for the RGI operator with the polynomial fit. 
Filled symbols represent interpolated values at p = 2.6, 2.9 and 3.2. A comparison is made between the 
ordinary renormalization condition (circle) and the improved condition (triangle). 
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FIG. 20: Scaling behaviors of RGI quark masses and continuum extrapolations, u, d quark mass and strange 
quark mass from the K input are defined by m q + m ros . 
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FIG. 21: Scaling behaviors of renormalized quark masses in MS scheme and continuum extrapolations. 
u, d quark mass (upper panel) and strange quark masses from the K input (lower left panel) and the <f> 
input (lower right panel). Results of this paper are represented by filled circles, together with the CP- 
PACS result with the perturbative renormalizationQ (solid squares). The VWI(AWI) quark mass with the 
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FIG. 22: Constraint bands in the CKM triangle from €k- The solid lines are central value and one standard 
deviation of the constraint with our Bk ■ The dashed lines are results with Bk adopted by the CKM fitter 
group [42] • For other inputs we use those given by the CKM fitter. 
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FIG. 23: Effective masses of the pseudo-scalar meson as a function of the temporal distance t at each quark 
mass. Lines represent the central value of the exponential fit of the propagator and its fitting range. 
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FIG. 25: Pscudoscalar meson mass squared (left) and vector meson mass (right) as a function of bare quark 
mass mfd. Lines show linear fits. 
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FIG. 26: Ratios of the matrix element to the vacuum saturation (|I.2|) as a function of t. 
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